Thermodynamics, spectral distribution and the nature of dark energy 
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Recent astronomical observations suggest that the bulk of energy in the Universe is repulsive and 
appears like a dark component with negative pressure (cj = Px/ Px < 0). In this work we investigate 
thermodynamic and statistical properties of such a component. It is found that its energy and 
temperature grow during the evolution of the Universe since work is done on the system. Under the 
hypothesis of a null chemical potential, the case of phantom energy {u < —1) seems to be physically 
meaningless because its entropy is negative. It is also proved that the wavelengths of the ti;-quanta 
decrease in an expanding Universe. This unexpected behavior explains how their energy may be 
continuously stored in the course of expansion. The spectrum and the associated Wien-type law 
favors a fermionic nature with u naturally restricted to the interval —1 ^ w < —1/2. Our analysis 
also implies that the ultimate fate of the Universe may be considerably modified. If a dark energy 
dominated Universe expands forever, it will become increasingly hot. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



The so-called dark energy or quintessence is believed to 
be the first observational evidence for new physics beyond 
the domain of the Standard Model of Particle Physics. Its 
presence, inferred from an impressive convergence of ob- 
servational results along with some apparently successful 
theoretical predictions, not only explains the current cos- 
mic acceleration but also provides the remaining piece of 
information connecting the inflationary flatness predic- 
tion with astronomical data Q . 

In spite of its fundamental importance for an actual 
understanding of the evolution of the Universe, the na- 
ture of this unknown energy component (which seemingly 
cannot be unveiled from background tests ) constitutes 
one of the greatest mysteries of modern Cosmology and 
nothing but the fact that it has a negative pressure (and 
that its energy density is of the order of the critical den- 
sity, ~ 10~^^ g/cm'^) is known thus far. This current 
state of affairs brings naturally to light some important 
questions (among others) to be answered in the context 
of this new Cosmology: 

1. What is the thermodynamic behavior of the dark 
energy in an adiabatic expanding Universe? Or, 
more precisely, what is its temperature law? 

2. How do dark energy modes evolve in the course of 
the cosmic expansion? 

3. What is the dark energy frequency spectrum? 

4. What is the thermodynamic fate of a dark-energy- 
dominated universe? 

In order to answer the above questions, it is neces- 
sary to make some fundamental hypothesis concerning 
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the unknown nature of the dark energy. In principle, the 
thermodynamic behavior and the fate of cosmic evolution 
are somewhat entertained and depend fundamentally on 
the previous knowledge of the intrinsic nature of the ex- 
panding components. 

In this letter we investigate the thermodynamic and 
some statistical properties of the dark energy assum- 
ing that its constituents are massless quanta (bosonic or 
fermionic) obeying a constant equation-of-state parame- 
ter, uj = Px/px, where Px is the pressure and px its energy 
density. To be more precise, most of the thermodynamic 
results remains true even if the dark energy medium is 
formed by massive particles. 

In what follows, answers to the above questions are 
sought by using only the observational piece of informa- 
tion that for dark energy lu is a negative quantity. Since 
Lu can also be seen as a continuous parameter, the case 
Lu — 1/3 (blackbody radiation or neutrinos) provides a 
natural check for our general results. Throughout this 
paper it will be assumed that the Universe is described 
by the homogeneous and isotropic Friedmann- Robertson- 
Walker (FRW) geometry 

ds^ ^^dt^~R\t){^^+r^de^+r^sin^ed(f'), (1) 

where k = 0, ±1 is the curvature parameter and R{t) is 
the cosmological scale factor. 

1. Dark Energy and Thermodynamics. To begin dis- 
cussing the first question, we recall that the thermody- 
namic states of a relativistic simple fiuid are character- 
ized by an energy momentum tensor T"'^ , a particle cur- 
rent iV" and an entropy current S"". By assuming that 
the dark energy is a perfect fluid such quantities are de- 
flned and constrained by the following relations Q 

= T"^;/3 = 0, (2) 
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Ar" = nw", iV";„=0, (3) 



S°' = nau°', 5'";a=0, (4) 

where (;) means covariant derivative, n is the particle 
number density, a is the specific entropy (per particle) 
and the quantities Px, Pxi n and a are related to the 
temperature T by the Gibbs law 

nTda — dpx — ^-^—^dn. (5) 
n 

By taking n and T as independent thermodynamic vari- 
ables and by using the fact that da is an exact differential, 
it is straightforward to show that the temperature evo- 
lution law is given by (IJ (a dot means comoving time 
derivative) 

or, equivalently, for ^ 

n = constr- ^ T^/^v = const., (7) 

since n scales with V"^, where V is the volume of the 
considered portion within the fluid. In the case of black- 
body radiation {uj — 1/3), the above equations yield the 
well known results n cc and T^V = const. It means 
that blackbody radiation cools if it expands adiabati- 
cally, a typical behavior for fluids with positive pressure 
{lo > 0). Physically, this happens because each portion 
of the fluid is doing thermodynamic work at the expenses 
of its internal energy. ^From Eq. Q, however, we find 
that dark energy becomes hotter in the course of the cos- 
mological adiabatic expansion since its equation-of-state 
parameter is a negative quantity. A physical explanation 
for this behavior is that thermodynamic work is being 
done on the system (see Fig. 1). In particular, for the 
vacuum state (w = —1) we obtain T oc V. Therefore, 
any kind of dark energy with negative pressure (Includ- 
ing the so-called phantom energy, u; < —1 be- 
comes hotter if it undergoes an adiabatic expansion. Nat- 
urally, this unexpected result needs to be explained from 
a more fundamental (microscopic) viewpoint, which is 
closely related to the remaining questions. However, be- 
fore discussing these questions, it is worth noticing that 
by combining Eqs. ^ and Q one obtains the general- 
ized Stefan-Boltzmann law (see for some independent 
derivations) 

Px{T) ^ V^t"^ , (8) 

where rjuj is a w-dependent constant. For lo = —1 one 
finds px = constant, as expected for the cosmological 
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FIG. 1: Dark Energy and thermodynamic work. The total 
energy and temperature of the dark component grow in an 
adiabatic expansion because work has been done on the sys- 
tem (_Ej; = p:cV cuT). From Eqs. (7), (8) and (9) we see that 
a vacuum filled Universe (a; = —1) is a limiting case where the 
temperature increases, the energy density remains constant, 
and the entropy is zero. 

constant case. The entropy of this ti;-fluid is also of in- 
terest. If the chemical potential is null (as occurs for 
uj — 1/3), the expression a = Sx/N — [px +Px)/nT de- 
flnes its specific entropy Therefore, the dark energy 
entropy can be expressed as 

Sx{T,V)^r^^{l+Lo)T^'''V, (9) 

from which the adiabatic temperature evolution law |(7J) is 
readily recovered. The above result completes our ther- 
modynamic description for the dark component. Note 
that the vacuum entropy is zero {uj = —1) whereas 
for phantom or supernegative dark energy (lu < — 1) 
the entropy assumes negative values being, therefore, 
meaningless. If the hypothesis of a null chemical po- 
tential is reasonable, this latter result poses a serious 
problem to a phantom energy description with basis on 
the usual XCDM parameterization and introduces a new 
thermodynamic lower limit to the dark energy equation- 
of-state parameter, i.e., to > —1. Such a limit is in fully 
agreement with theoretical arguments that motivate the 
so-called A barrier and goes against several observa- 
tional analyses which seem to favor phantom cosmologies 
over ACDM or usual quintessence scenarios 0, H, 01- It 
should be recalled that thermodynamic states with nega- 
tive pressure are metastable and usually connected with 
phase transitions (the most known exam ple at the level 
lab is an overheated Van der Waals liquid [i3)- It should 
be stressed, however, that the thermodynamic behav- 
ior derived here holds regardless of the specific physical 
mechanism responsible for the negative pressure. 
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2. Effects of the Expansion on the Dark Energy Modes. 
In what follows, we hipothcsize that the constituents 
of the dark energy are massless particles (bosons or 
fermions) and that the differences between this fluid and 
blackbody radiation (or neutrinos) are macroscopically 
quantified by the equation-of-state parameter lo. As we 
shall see, this hypothesis leads to a microscopic picture 
in harmony with the above thermodynamic behavior. 

The adiabatic theorem guarantees that if a hol- 
low cavity (e.g., our Universe) containing dark energy 
changes adiabatically its volume, the ratio between the 
energy of a given mode and the corresponding frequency 
remains constant, that is, E^/v = const., for any proper 
oscillation. Thus, if px(T,v) is the spectral energy den- 
sity inside an enclosure with volume V at temperature T, 
this adiabatic invariant guarantees that Px{T^ i')dvV/v — 
const. Now, recalling that the energy density in the band 
dv varies with the temperature in the same manner as 
the total energy density, one may write [see ijSJl] 

T^y/z/ = const. (10) 

and, since Tiv = const., it follows that T/v is an in- 
variant. Therefore, whether a hollow cavity containing 
dark energy is expanding adiabatically, the wavelength 
of each mode satisfies 
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FIG. 2: Pressure effect on the wavelength at two different 
times (ti, t2) in an expanding Universe. The wavelengths from 
radiation with positive (negative) pressure increase (decrease) 
in the course of the expansion. Energy is continuously stored 
in the dark energy modes (p < 0). This unexpected behav- 
ior is a consequence of Eqs. O and IIH . It explains why 
the internal energy of the dark radiation increases during the 
evolution of the Universe. 



AT = const. (11) 

The above results allow us to understand from a micro- 
scopic viewpoint why the net energy of a dark component 
increases during the expansion. As shown in Figure 2, the 
wavelengths of any radiative fluid with positive pressure 
{uj > 0) increase in virtue of the expansion thereby low- 
ering the total energy in accordance with the generalized 
Stefan-Boltzmann law. However, if w < 0, we have seen 
that the temperature grows whether the fluid undergoes 
an adiabatic expansion [see Eq.Q], and this is accom- 
panied by a decreasing in each wavelength A in order to 
satisfy the above relation. This, therefore, constitutes a 
microscopic explanation of why the energy in a given por- 
tion of the dark component increases (see Figure 1). In 
particular, the vacuum state, i.e., uj — —1, behaves as a 
limiting case for which the total energy increases during 
the expansion while its energy density remains constant. 

3. The spectrum of dark energy. Let us now quantify 
the influence of a negative pressure on the general form 
of the spectral distribution. This is an important point 
because relations ^ and JSJ must be recovered from the 
frequency spectrum. To that end, we consider again an 
enclosure containing dark energy at temperature Ti and 
focus our attention on the band AAi centered on the 
wavelength Ai whose energy density is Px{Ti, Ai)AAi. If 
the temperature Ti changes to T2 due to an adiabatic ex- 
pansion, the energy of the band changes to Px{T2, A2) AA2 
and, according to Ea. (|ll|l . AAi and AA2 are related by 
AA2/AA1 = T1/T2. Now, since one can assume that dis- 



tinct bands do not interact, it follows that 

P:,(r2,A2)AA2 ^ ,T2 ^ 

p,(Ti,Ai)AAi " ■ ^ ' 

By combining the above results, and using again the adi- 
abatic invariant (|ll|l . we obtain for an arbitrary com- 
ponent px{T,X) — a A ~4>{XT), or still, in terms of 
frequency 

Px{T,iy) = aiy^(j}{^), (13) 

where a is a constant with dimension of 
[Energy. Lenght^'^. Time ''"'^^'^''^'^], and (j) is an arbi- 
trary dimensionless function of its argument. The above 
expression is the generalized Wien-type spectrum for 
dark energy. 

The specific form of the dimensionless function <j){v /T) 
can be determined using different approaches, among 
them, the arguments put forward by Einstein in his 
original deduction of the blackbody distribution [T^ . 
For a massless bosonic gas, for instance, (f>{h'/T) is ex- 
actly the occupation number (without chemical poten- 
tial), and the spectrum takes the final form px{T^v) = 
av^ X \e^^/^BT _ ^ggg fpj. uiore detailed discus- 
sion). Therefore, in order to include the case of fermions, 
we write the spectral distribution (13) as 
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which is the most natural extension of a Planck (or Fermi- 
Dirac) spectrum describing the a;-family of dark energy 
radiation. Einstein's derivation follows for w = 1/3 [l4| 
and, more important, the macroscopic relations O and 
are recovered from the above spectral distribution, 
i.e., 



niT)= r P^^dj^^n^T^, (15) 
Jo hv 



and 



Px{T) = / Px{T, v)dv = r]^ T~ 
Jo 



(16) 



where the pair of constants {fj^,rii^) depend on a, ks 
and uj, and also on the underlying statistics obeyed by 
the dark quanta. 

Another interesting point is related to the wavelength 
for which the distribution (|14(l attains its maximum 
value. It is determined by the algebraic condition 



e"'' ± 1 = ±uj{l + 2uj)-^x, 



(17) 



where x = hc/ksXT and the signs ± stand for the 
fermionic and bosonic cases, respectively. The displace- 
ment Wien's law now reads 



he 



1.438 



kBx'{Lj) x'{lo)' 



(18) 



where x'{u!) is the root of the above transcendent equa- 
tion. As one may check, in the case of bosons, a large 
set of nontrivial solutions exist with positive tu, but all of 
them are out of the dark branch (w < 0). However, by 
assuming that the basic constituents of dark energy are 
massless fermions, the existence of solutions fixes the up- 
per limit Lu < —0.5. In particular, for uj — —2/3, a value 
compatible with many astrophysical data P, , we find 
X^T = 1.943 cm.K while for uj = -1, A„T = 1.123 
cm.K. 

4- Dark energy and the thermodynamie fate of the ex- 
panding Universe. The fate of our Universe, i.e., whether 
it will eventually re-collapse and end with a Big Crunch, 
or will expand forever, is a matter of great interest among 
cosmologists (see, e.g., ^^). In the standard view, if 
the Universe expands forever it will become increasingly 
empty and cold. However, whether our analysis provides 
a realistic description of the the dominant dark energy 
component, such a destiny is not so neat because of its 
strange thermodynamic behavior. 

In the context of a FRW-type geometry {V oc R^), Eq. 
{71) implies that T cc R~^'^ or, equivalently, that the gen- 
eral redshift-temperature law reads Tx{z) = T°(l -I- z)^'^, 
where T° is the present value of the dark energy tem- 
perature. As discussed earlier, Tx{z) grows in the course 



of the expansion while the observed average temperature 
of the Universe is given by the decreasing temperature 
of the CMB photons. Therefore, an important point for 
the thermodynamic fate of the Universe is to know how 
long the dark energy temperature will take to become the 
dominant temperature of the Universe. A basic difficulty, 
however, is that the present-day dark energy temperature 
has not been measured and, moreover, the a parameter 
appearing in the spectrum [Eq. (14)] (or, equivalently, 
the r]c^ in the expression of the energy density) is still 
unknown [l7| | . A very naive estimate of T° can be done 
by considering the observational evidence for the density 
parameters, fix ^ 0.7 and ^Ir ^ 10~^ (radiation). For 
UJ ^ —1 one finds T° ~ (10^a/77i^)T+^ , where a is the ra- 
diation constant. Therefore, any estimate of the present 
dark energy temperature depends crucially on the ratio 
parameter a/rj^. In particular, for co = —2/3 we find 
T° ^ 10^^"(77_2/3/a)^, so that temperatures of the order 
of 10~^K are easily obtained for 77-2/3/0 ^ 10^. Albeit 
many aeons might be necessary for the Universe to enter 
in the dark energy temperature regime, the only possi- 
ble conclusion is that the ultimate fate of the Universe 
may be considerably modified: a dark energy dominated 
Universe expanding forever will become increasingly hot. 

Summarizing, we have proposed that the constituents 
of dark energy are massless particles (bosons or fermions) 
whose collective behavior resembles a kind of radiation 
fiuid with negative pressure. Through a thermodynamic 
analysis and basic statistical considerations we derived 
some relevant properties obeyed by this mysterious com- 
ponent. As expected for a consistent treatment, all ther- 
modynamic relations are recovered from the statistical 
approach. Perhaps, more important, with basis on a 
Wien-type law it was possible to show that a fermionic 
nature is clearly favored because there are nontrivial so- 
lutions within the dark energy branch {lo < 0). It was 
also proved that the existence of such solutions requires 

to be < —0.5. When combined with the limit obtained 
from the entropy calculation [Eq. (9)], this result con- 
strains the equation-of-state parameter to the interval 
— 1 < UJ < —0.5, which is surprisingly close to the range 
required from many different astrophysical data P, 0| . 
Such results means that some kind of massless particles 
obeying the Fcrmi-Dirac statistics should be considered 
among the candidates for dark energy proposed in the 
literature. However, it should be stressed that our treat- 
ment does not eliminate the case of massive fermions for 
which most of the thermodynamic properties discussed 
here remains true. 
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